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Abstract 

We analyze the zero energy solutions, of a two dimensional system 
which undergoes a non-radial symmetric, complex potential V {r,(j)). By 
virtue of the coherent states concept, the localized states are constructed, 
and the consequences of the imaginary part of the potential are found 
both analytically and schematically. 



1 Introduction 

Since the early years of quantum mechanics the exact solvability of quantum me- 
chanical models have attracted much attention. Some exactly solvable models 
have already become typical standard examples in quantum mechanical text- 
books. However, it was believed that the reality of the spectra of the Hamil- 
tonians, describing quantum mechanical models, is necessarily attributed to 
their Hermiticity. It was the non-Hcrmitian PT-synimetric Hamiltonians pro- 
posed by Bender and Boettcher [1] that relaxed the Hermiticity condition as 
a necessity for the reality of the spectrum [1-7]. Herein, P denotes the par- 
ity (PxP — —x) and the anti-linear operator T mimics the time reflection 
(TiT = Recently, Mostafazadeh [8] has introduced a broader class of non- 
Hermitian pseudo-Hermitian Hamiltonians (a generalization of PT-symmetric, 
therefore). In these settings [8-19], a Hamiltonian H is pseudo-Hermitian if 
it obeys the similarity transformation: j]Hj]^^ — where is a Hermitian 
invertible linear operator. On the other hand, the study of the E = bound 
states have found many applications in various fields [20-24]. Long ago, Daboul 
and Nieto [25,26] have discovered that, an attractive radial power low potential, 
V{r) ~ r~'^ for v < —2 and v > 2, passes through the E = normalizable 
solutions. More recently Makowski and Gorska, established the classical cor- 
respondence localized states of a system with zero energy and a general form 
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of power low potentials [27]. In their work, it was shown that the classical 
trajectories of the particle precisely matched with the localized quantum states. 

In this work we attempt to understand the consequences of adding an imag- 
inary term to a potential whose zero energy level passes through the E = nor- 
malizable bound state solution (namely V{r) = —Tr"^). We believe that this 
kind of study is necessary, to find relations, if any between the non-Hermitian 
quantum mechanics and the classical mechanics. 

The organization of our paper is as follows: In Sec. (2) we give an analytical 
solution to the Schrodinger equation of a zero energy particle under our cho- 
sen complex potential. We continue in Sec. (3) by adapting a closed form of 
the localized states from literature and then we set up and plot the classically 
equivalent coherent states of the system. We conclude our paper with Sec. (4). 



2 Analytic solution of the Schrodinger equation 

Two dimensional Schrodinger equation for a zero energy particle under a com- 
plex, non-radial symmetric potential is given by 
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V'(r,</.)=0 (1) 



where 

y^rA) = ~-^,e^'' (2) 

and for our purpose, F and A arc some non-negative real constants. Before we 
go further, finding out the symmetric properties of V (r, (j)) and consequently the 
Hamiltonian of the system may give some connections between this potential 
and the well known 'PT-symmetric type potentials which are studied in the 
literature. 

Let us introduce an operator 6 which is defined as 

e : i ^ -1,0^ 27r- (3) 

in which i = y/^-i and (j) is the usual azimuthal angle. One can easily show that 
is a non-Hermitian, invertible operator whose inverse is given by 

e-i = e (4) 

and it can be decomposed into the two other operators T and 11 such that 

e = nr (5) 

in which the definition of these operators are given by 

n : </> ^ 27r - (/) (6) 
T : -i. 
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n is Hcrmitian and invcrtiblc such that 11""'^ = 11 = H'^, and T is the usual 
time reversal operator. It is remarkable to observe that the Hamiltonian of the 
particle under the potential (2) is IIT-symmetric which means 



H = H 



nr 



(7) 



Also it is said that the HT-symmetry of a Hamiltonian H is unbroken if all 
of the eigenfunctions of H are simultaneously eigenfunctions of UT. It is easy 
to show that if the IIT-symnietry of a Hamiltonian H is unbroken, then the 
spectrum of H is real (to see the proof, one may see [30]). 

We come back to the equation (1) and as usual, we take 

V'(r,,^) = i?(r) $(</,) (8) 
and after making substitution, in Eq. (1) we get a set of two equations as 



ld_ 
r Or 



d_ 
dr 



— + 

and 



2mT 



R{r)=0 



ft' 



-/2$(<^) 



where I is a constant to be identified. 

By introducing a new dimensionless variable p as 



P = 



(9) 



(10) 



(11) 



where ao is a positive constant, the Eqs. (9) and (10) become 



in which 



1 d 



pdp \ dp 
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and 



R{r)=0 



+ A^e'*$ {(I)) = -V^ {(/)) . 



7 = 



2mT 

and 
2mA 



(12) 



(13) 



(14) 



(15) 



In the angular part of the Schrodinger equation we change the variable, and 
introduce 

X = (16) 
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hence the Eq. (13) reads 

This is the modified Bessel ODE, such that its complete solution is well known 
as 

$(X) = Cil2i{2X^) + C2K2i{2X^) 

or (18) 
$ (0) = Cihi (2Ae*'^/2) + C2K21 (2Ae^'^/2) 

in which (z) and -ftT^ (z) are the modified Bessel functions. 
This solution should satisfy the following boundary condition 

^ {(/,) = ^ {(I, + 27r) (19) 

or equivalently 

Cihi (2Ae'^/2) + ^2X2; (2Ae^'>/2^ = (20) 
Chi (2Ae^(^+2^)/2) + C2i^2; {2Xe'^f+^^^/') . 

Since I,^ (z) and K,^ (z) arc two independent solutions of the modified Bessel 
ODE , a class of solution is possible when we put C2 = 0. Therefore we get 

I21 (2Ae^^/2) = I21 (2Ae^(^+2")/2) . (21) 

In accordance with the following property of the modified Bessel functions 

[28] 

(ze"*") = e"""'"'!^ {z) (22) 

where ^; is a complex variable, m is an integer and v is a, real number, one can 
choose 

z = 2Ae*'^/2 (23) 

and m = 1 in Eq. (22) to get 

I21 {ze}^) = e'^^^hi (z) . (24) 

This equality is valid for all z and I in their domains, but if one considers / to 
be an integer, we will get 

I21 {ze:^) = hi {z) (25) 
which is equivalent to Eq. (22). Therefore / is found to be an integer i.e. 

Z = 0,±1,±2,... (26) 
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Of course as a different possibility, one can choose Ci = to find a different 
class of solution but the following property [28] 

(ze'"''') = e-"''"''K^ (z) - m sin mvK esc vt^U {z) (27) 

it comes to our setting as 

K.21 {ze") = e-''^^'K2i {z) - wi sin 2ln esc pwhi (z) (28) 

which obviously does not admit any solution. 

As a result, the solutions of the angular part of the Schrodinger equation 
can be written explicitly as 

= QA/2i (2Ae*"^/2) (29) 
/ = 0,1,2,... 

where, since /_„ (z) = In {z) , we just consider I to be non negative, and Ci\ are 
the normalization constants given by 



One may note that I still can be interpreted as the angular quantum number, 
since 

< L>=<$i|-i/i|-|$,>= (31) 



< 



Cahi (2Ae'^/2) | - ^fi^lQxhi (2Ae^*/2^ >= 



>= ih. 



The radial part of the Schrodinger equation, can be considered as the Bessel 
ODE if one defines 

£ = i (32) 

and therefore the Eq. (12) becomes 

f ^i?(0 +^1^(0 + (7¥ - l')R{0 = (33) 
which admits a complete solution 

i?i^(0=CiJ;(70+C2Vi(70- (34) 

A physical, normalizable solution which for Z > 1 corresponds to the bound 
state is given by [25,26] 

RiAC) = Ni^JihO (35) 
/ = 2,3,... 
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in which Ni^ are normahzation constants given by 



_ 2 (1 + 1)! 

Now we are ready to write the complete solution of the Schrodinger equation 
(i.e. wave function), by using the Eqs. (29) and (35) as 



= a.^^f±|/.(2A.*=).,(2^.). ,37, 

We notice that, with ^ > 1, the only complex part of ij^i ^ ■y{r^(j)) is the 
modified Bessel function. As a matter of fact, the effect of 8 (introduced in Eq. 
(3)) on ^/j; ^ (r, (f) is equivalent to the effect of on I21 (2Ae*'^/^) . Therefore, 
by using the expansion form of the modified Bessel function, this equation can 
be written as 



■2(s+l) 



in which, clearly this is invariant under 9 (i.e., V';,7,a ^) = V';,7,a (^j 'P))- 

2.1 A realized approach to the problem 

In this section we will consider the following radial symmetric real potentials 

^±W = -^±^ (39) 

where V and A arc some positive constants as before. One should notice that 
the negative branch of the above potential is same as the potential in Eq. (2) in 
an attractive form, and the positive branch of it is same but in a repulsive form. 
The Schrodinger equation (1). after the usual separation method and change of 
variable, with the potential (38) comes to a set of two separated equations as 



pdpYdp) p^^ p^ 



R{p)=0 (40) 



where 



and 

= (41) 

d(p 



in which the positive (negative) sign is related to the +^ ( — ^),and the other 
factors are defined as before. One can easily show that, the final solution of the 
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Schrodinger equation with the potentials (38) can be written as 



in which / must be greater then 1 to have bound states. In the sequel we will 
use the closed forms of the infinite number of degenerate wave functions (i.e., 
these states have same energy equal to zero), presented in the Eq.s (37) and 
(41) to construct the classical correspondence localized states. 



; + l ! 



(42) 



3 Zero energy localized states 

In references [25-27] it was shown that the trajectory of a classical particle which 
experiences a real potential in the form of the F— part of the potential considered 
in this work (2) (i.e. — ^) is given by 

a; = |(l + cos((/.-(/.J) 

y = § sin (</, - 0J ^^-"l 

which represents a circle with radius a/2 (i.e., if one set (j)^ = 0, this becomes 
{x - a/2f = {a/2f) such that 

a = ^J^mFfl? (44) 

and L is the conserved angular momentum of the particle. We notice that, the 
classical correspondence localized state (see [27] and the references therein) of 
the potential (2) while A — > must have probability peak in accordance with 
the classical trajectory (i.e., a circle with radius a/2 as implied by Eq. (42) [27] 
). Our aim in the following is to see the effect of the A— part of the potential (2) 
on the shape of the classical correspondence localized states. To this end, first 
we find the localized states of the original potential (2) in a closed analytical 
form, then we follow similarly but for the case when the potential is in the forms 
ofEq. (38). 

An available method to derive the corresponding localized states by using the 
solutions of the Schrodinger equation given in the previous sections, is based on 
the concept of deformed oscillator algebras [27,29]. Therefore a suitable explicit 
form of the localized state over the infinite number of degenerate states ( with 
E = 0) tpijx (^) 0) reads (see [27,29] and the references therein) 

^ /7V\ 1/2 



in which k = 1 — 2 and r = Ae*^° , where A and 9o are some real constants. 
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3.1 Results 



The classical correspondence localized states of a particle undergoes the poten- 
tials (2) and (38), by choosing N = 7 may be written as 

8^ E C'(2+fe)A\/^^%^-f(4+2fc) hXe'-) J(2+fc) (i) V = Vc 
V fe=0 V V / 

(46) 

in which and A are set to be one and for convenience Vc and V± refer to the 

potentials (2) and (38) respectively. Within the Fig.s (1) to (7) some density 
plot of \^^\ are given in terms of different values of A. In the Fig. (1) we plot 
with A = (i.e., Vc = V± = —-pc) and the classical trajectory of the 
particle ( this figure was reported in the Ref.[27]). In the Fig.s (2) to (4), (a), 
(b) and (c) refer to the potentials V-, V+ and Vc, respectively. In the Fig.s (5) 
and (6), (a) refers to V+ and (b) refers to Vc- Finally Fig. (7) refers to 14 . 



3.2 Behavior of the complexified modified Bessel func- 
tions 

In this section we give a short description of the behavior of the complexified 
modified Bessel functions hiiz) to explain why with V = Vc the plot of the 
probability density with large values of A are very localized around (/> = 0. It is 
not difficult to see that, in the case oiV = Vc the A— part of the original potential 
(2) goes into the (p-pa-it of the Schrodinger equation and therefore the entire 
effect of this term appears in the 0-part of the wave function. Therefore the 
contribution of Cixhi (2Ae"^/^) in the final form of the wave function, instead 
of the usual 0— part in the wave functions (LCj-^t^e*'*^) of the case of y = V± 
is the reason of the great localization about small 0. Let us write 



(2Ae^^/^) 



/2i 2Acos^ +i2Asin^ (47) 



which shows that the square root of the real (imaginary) part of the potential 

directly goes through the real (imaginary) part of the argument of I21 (2Ae*'*/^) . 
To see the behavior of I21 (2Ae*'^/^) in terms of the real (imaginary) part of its 
argument we rewrite the last equation as 



i2l 



hi 2Ai cos ^ + i2X2 sin ^ (48) 



which in the limit of Ai = A2 = A turns out to be hi (2Ae*'^/^) . Figures (8), (9) 



and (10) show that once A2 vanishes \hif does not change much, but once Ai 

becomes zero, \hi\^ decreases strongly. Also once Ai takes a larger value (it 
does not matter what is the value of A2), \hi\'^ takes much higher value close to 
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(f) ~ or Itx. Wc conclude therefore that, the imaginary part of the argument of 
1-2X (and imaginary part of the potential therefore) does not contribute much in 
the localization of the around small </>. In contrast, the real part, and its 
0— dependent part (i.e., cost/)) causes such a great localization. 

4 Conclusion 

In conclusion, we concentrate ourselves and wish to comment on the figures. Fig. 
(1) is our reference figure, i.e., the classical correspondence localized state when 
the A— part of the potential vanishes. Fig.s (2-4) clearly show that the effects 
of the A— part in the Y^. and Vj^ cause a higher localization while in the we 
see a lower localization. It is remarkable to observe that, the magnitude of the 
effects (within these three cases) Vc is greater than others. Fig.s (5,6) show that 
as A takes large value the radius of the localized state corresponding to the V+ 
decreases while the distribution of the |*7|^ does not change. For the case of 
potential Vc the radii of the localized states are fixed while the distribution of 
the probability density \^-!^ is changed so that the particle seems to be localized 
more around = 0. In the Fig. (7) we see the effect of the A— part in the Vc as 
a great localization about = 0. Obviously the figures imply that in the case of 
V = Vc, the particle is localized around = 0, which is a direct consequence of 
the A— part of the potential. We notice that, our approach to the problem is in 
a closed analytical form, where all numerical results are based on the analytical 
solutions. There is no need to comment that any other approach will definitely 
lack the advantages of an exact analytical solution. 
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Figure caption 

Figure (1): A plot of probabihty density \^ n (r, 4>)^ for = 7 (i.e. 2 < I < 
9), 6*0 = 0, A = 1, e = (i.e. the V{r, (p) is real) and A = 0.0. This is the localized 
state corresponding to the classical trajectory of a particle which experiences 
just the first term of the potential, i.e. V{r,4>) = —T/r^ and therefore it is the 
reference plot. 

Figure (2): A plot of probability density |*jv (r,0)|^ tov N = 7 (i.e. 2 < 
I < 9), 9o = 0, A = 1, X = 0.1. Also (a), (b) and (c) are correspondence with 
V{r,(j)) = V-,V+ and Vc respectively. 

Figure (3): A plot of probability density |*jv {r,(p)f N = 7 (i.e. 2 < 
I < 9), 9o ^ 0, A = 1, \ = 0.5. Also (a), (b) and (c) are correspondence with 
V{r, (j)) = V-,V+ and Vc respectively. 

Figure (4): A plot of probability density |\l/Ar (r, (/))|^ for N = 7 (i.e. 2 < 
I < 9), 9o = 0, A — 1, X = 1. Also (a), (b) and (c) are correspondence with 
V{r,(p) = V-,V+ and Vc respectively. 

Figure (5): A plot of probability density [^at (r, (/))|^ for AT = 7 (i.e. 2<l< 
9), 0o = 0, A = 1, A = 5. Also (a) and (b) are correspondence with V{r, (p) = V+ 
and Vc respectively. 

Figure (6): A plot of probability density [^at {r,(j))\ for N = 7 (i.e. 2 < 
I < 9), 6o = 0, A = 1, X = 10. Also (a) and (b) are correspondence with 
V{r, (p) — y+ and Vc respectively. 

Figure (7): A plot of probabihty density [^-at (r, (/))|^ for A/' = 7 (i.e. 2<l< 
9), 00 = 0, A=l, X= 100. 

Figure (8): A plot of hi (2X1 cos (^|) + i2X2 sin j j in terms of (p, for 
some different values of Ai and A2. 
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Figure (9): A plot of hi (2A1 cos (|) +i2A2sin (^|) 
some different values of Ai and A2. 

Figure (10): A plot of hi (2X1 cos (|) + i2A2 sin 
some different values of Ai and A2. 



in terms of 6, for 



in terms of 6, for 
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